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Abstract 

We study the generation of entanglement for interacting cold atoms in an optical lattice. The 
entanglement is generated by managing the interaction between two distinct atomic species. It 
is found that the current of one of the species can be used as a good indicator of entanglement 
generation. The thermalization process between the species is also shown to be closely related to 
the evolution of the current. 
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Entanglement is an important ingredient for the development of quantum information 
processing. Its study has been stimulated by its potential implementation in quantum 
computing aigorithms, making it ubiquitous in quantum systems, becoming particu.arl y 
interesting in many-body systems Fast information processing is expected to involve 
many elements in the system, making particularly attractive the potential use of ultracold 
atoms. On the other hand, cold atom systems have proven to be ideal devices for ultra- 
high precision measurements j^j]. Along these lines, several works involving entanglement in 
ultracold atoms have been put forward ^-^]. 

In the present work we study a system where entanglement not only can be efficiently 
generated, but also accurately estimated using current of particles. Here we consider that 
the parts of the system are two sets of different species of atoms moving on a lattice with 
tunable interaction. A sudden interaction quench between the two distinct atomic species is 
applied in order to generate entanglement. In cold atom systems interaction can be easily 
manipulated, representing a great advantage with respect to semiconductor superlattices. 
Such a freedom to manage the interaction has been demonstrated experimentally |6| and 
has proven to be fundamental in phenomena such as coherent destruction of tunneling in 
a many-body system ?]], control of atomic localization along a lattice jsj, stabilization of 
Bloch oscillations for a train of solitons in Bose-Einstein condensates (BEC) [9|, among 
others. On the other hand, increasing the interaction between atoms modifies the flux of 
particles along the lattice, giving rise in the limit of large interaction to a quantum phase 
transition from a superfluid state to a Mott insulator [10] . Studies involving atoms moving 
along a ring shaped optical lattice have exploited such interaction combined with magnetic 
fields to generate a current of particles 111] , as well as for the implementation of an AC 
quantum motor jl2 |. 

We study a system where initially two sets of different species move independently along 
a ring with a finite current. As the interaction is turned on it not only modifies the currents, 
but also induces entanglement between both groups of atoms. We investigate the correlation 
between current and entanglement where the detection of the particle current variation 
provides a good measure of entanglement generation. Interacting species of condensed atoms 
can be achieved either by using two kind of atoms, such as 87 Rb and 40 K [13] or from the 
same isotope with different internal states [5|. The Hamiltonian for two interacting species 
of ultra-cold atoms reads 
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H = H A ®i B + H B ®l A + H int , (1) 
where the Hamiltonian for each specie D = A, B is described by the Bose-Hubbard model 



3=1 3=1 

with dj (J = a, 6) the creation and annihilation operator for the particles at site j, = d^dj 
is the particle number operator and L is the total number of sites. The first term describes 
the hopping or tunneling between adjacent sites in the lattice with a tunneling strength C. 
Here the hopping appears modified by a twist factor e % ^l L which is generated via a magnetic 
field and produces a finite current along the ring 10[ . The above tigh-binding Hamiltonian 



is valid within the weak-coupling regime, that is, for optical lattices with a depth > 5E. 



where E r = h 2 k 2 /2m is the single photon recoil energy and m is the atomic mass 14]. 
Hereon we assume the same mass for both atomic species. The second term characterizes 
the on-site interaction with strength V which, for simplicity, we consider equal to zero in 
what follows. 

The interspecies interaction Hamiltonian is given by 

L 

H int = Uj2ni®nf, (3) 

3=1 

where U denotes the interspecies interaction strength. Interspecies interaction can be ma- 



nipulated around Feschbach resonances [13J allowing great independence in tuning the inter- 
actions. Management of interactions provides a valuable tool for the control of the dynamics 
including generation of entanglement. Since we are interested in how the entanglement is 
generated between the species, we investigate the scenario with tunable interspecies interac- 
tion U(t). The full Hilbert space is spanned by the direct product of the single Fock states 
\ m A) <8> \^b) with dimension Na x iV^. For three lattice sites the dimension of the Hilbert 
space for one of the species is N D = {n D + l){n D + 2)/2, where n D is the corresponding 
number of atoms. 

Since we want to study the transport in a ring we use the periodic boundary condition 
\L + 1) = |1). The current operator for the particles in the ring is J = J a <8> 1b + 1^ <8> Jb, 
where Ja,b are the respective current operators for both species, defined as the sum of flux 



difference of particles between adjacent sites 



15| 



J ° = D^^+i - e-^ L d) +l d 3 ). (4) 
j'=i 

A good indicator of the degree of entanglement in a bi-partite pure state is provided by 
the Schmidt number /C defined as the reciprocal of the purity of the reduced density matrix 

H i 

Tr B{p 2 B ) 

For the sake of convenience, we define the normalized Schmidt number 

K = {Kq- 1)/(A - 1) where < K < 1 (6) 

and A = min(NA, Nb) is the dimension of the system, so the maximum degree of entangle- 
ment is found at /C = 1. 

We consider an L = 3 ring where the interspecies interaction is initially off and the cloud 
of particles is prepared in the ground state of the optical lattice by using adiabatic loading 
techniques. Particles are then set in motion by tuning a magnetic field [llj]. Under these 
circumstances the Hamiltonian is reduced to the tight-binding description where one can 
easily find the single particle eigenstates 

4i°) = ^E eimfe ^i°)' (?) 

v m=l 

where |0) is the vacuum state and k s = 2its/L, s — 0, 1,L — 1. Thus, the ground state of 
the non-interacting system with ha particles of type A and % particles of type B becomes 

I \^Jrnk A J ) 

3 (n A +n B )/2 ^ nA \ nB \ 



' \ / \ n B 

E^i (J2 ellkBb i) i°) 

k m J \ I J 



where {k B ) is the wavevector corresponding to the single particle ground state for A (B) 
particles. To get the two subsystems entangled, in analogy to experiments for the generation 
of entanglement using a nonlinear crystal, we consider the scenario where an attractive 
interaction strength is suddenly turned on at t = 0. Thus, we have a time-independent 
Hamiltonian before and after the interaction with constant interaction amplitude. 

Let us consider first the strong interaction limit \C/U\ — >■ 0, such that at times t > 
the dynamics becomes determined by the interspecies interaction H int . Thus, for a constant 
interaction, the wavefunction evolution is given by e~ t ~^ Lt \ijjo). From now on, we set a 
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small number of particles for the subsystem B, and an arbitrary number of particles for the 
subsystem A. For instance taking ub — 1 and n A = N results in the evolved wavefunction 

N 

\^ p imk B p ik A J2? =1 h p iL r(T,^Li ,m) 

3(iV+l)/2 x /jyT 



h,...,h,rn j=l 



Using Eq. Owe find the density matrix for the composite system p AB = \ip(t))(ip(t)\ which 
after tracing over the system A yields the reduced density matrix 



3 v+ 

nl+n2+n3=N ,m,m' 



D ik B (m-m') x 



e k ^ =l !j "" ^'^6^10) (0|6 W ,. (10) 

This result can then be substituted into Eqj5]to estimate the entanglement. Exact calcula- 
tions for the Schmidt number in the scenario of few particles can be found. In particular, 
for N = 1, we find 

27 

°~ 15 + 8cos(Ut/h)+Acos(2Ut/h) 1 ' 

which maximizes at 27r/3 + 2stt, and 47r/3 + 2sit with s = 0,1,2.... The evolution for 

the normalized Schmidt number is shown in Fig. Interestingly, a short time after the 

interaction is turned on the entanglement increases and reaches its maximal value. As shown 

in the figure, the entanglement remains high for a lapse of time before going down to zero 

again, repeating periodically with a time scale determined by U. If after some time the 

interaction is turned off, the entanglement of the system is preserved. Similar features are 

observed for = 2 in Fig. [T](b), where we also note an enhancement of the windows with 

high entanglement. 

Let us analyze now what happens with the current of the subsystem B at t > 0. Consid- 
ering again the evolved wavefunction \if>(t)) for A^ = 1 we find 

of 

Mt) = {J B ) = —sm(k B + 0) (1 + 2cos(Ut/H)) , (12) 
9 

which is a periodic function of the phase (p. For t > 0, the current becomes a time dependent 
function with a single frequency U. Interestingly, \Jb\ goes down as the entanglement 
increases (see Eq frTl) . Moreover, the maximal entanglement points exactly match with the 
zeroes of EqJT2l In fact, it is convenient to introduce the normalized current 



J(t) = 1 - 



J B (t) 

Mo) 



(13) 
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FIG. 1: (Color online) Evolution for attractive interaction U = — 10 4 C. Black line: Schmidt 
number 1C. Red line: Normalized current J . In both panels the exact analytical expression is 
superimposed, (a) one particle in A and one particle in B; (b) two particles in A and one particle 
in B. 



which shows surprising similarity with the evolution of entanglement JC(t) (see Fig. 1(a)). 
Such a result is interesting by itself, suggesting the use of current as a witness of the entan- 
glement. 

Contributions owing to the hopping elements in Eq. [1] have been neglected so far. To 
complement our previous analysis, the presence of a small contribution of the tigh-binding 
Hamiltonian is now considered within the strong interaction regime \C/U\ <C 1. As shown 
in Fig. Et^a), the entanglement evolution between one particle in each subsystem exhibits a 
similar behavior as for the zero hopping scenario, preserving the strong correlation between 
JC and J . Indeed, we observe again that the current vanishes at the points of maximal 
entanglement and that the absolute value of the current decreases as entanglement increases. 
However, increasing the number of particles enhances the role of the hopping term generating 
fluctuations of the previous profile, as depicted in Fig. E(b). Nevertheless, the normalized 
current of subsystem B still shows signs of strong correlation with the entanglement. 

We are interested in estimating the aforementioned correlation within a short timescale 
r s = 2tu and a much larger one T\ = 20tjj, where tjj = h/\U\ is the characteristic time of 
the system. To quantify the degree of correlation between the normalized current and the 
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FIG. 2: Left and right panels: (Color online) Upper graph: Normalized Schmidt number K, vs. 
time. Lower graph: J vs. time. The parameters are U = —IOC and (f> = n/10. Left: one particle 
in A and one particle in B. Right: two particles in A and one particle in B. 



normalized entanglement, we consider the Pearson product-momentum coefficient number 

R _ E[(lC-^)(J- f xj)] ^ 

where /i's and cr's are the mean values and standard deviations respectively. The absolute 
value of the coefficient R is close to unity when the variables involved are highly correlated 
and approaches zero when there is no correlation. Figl3](a) exhibits the results obtained 
from numerical simulations for two distinct sets of atoms number one and two particles in 
subsystem B and a finite hopping C = \U\/10. For ub — 1 the correlation between /C and 
J remains high even when the number of particles in subsystem A goes up to A^ = 17 for 
both time scales considered. Furthermore, for % = 2 a high degree of correlation appears 
to be almost independent of A" for the short time scale t s , whereas for a much larger time 
T\ the correlation decreases when the number of particles increases. 

Remarkably, for the scenario with large number of particles in A the entanglement does 
not decrease significantly after its initial growth as shown in Fig. [3](b) for AT = 11. Instead, 
the Schmidt number remains high, displaying only small fluctuations around its average 
value. This behavior arises because the finite hopping between adjacent sites lifts the degen- 
eracies of Hi n t such that the small subsystem B can now equilibrate with the much larger 
system A. In other words, particles of type A play the role of a thermal bath for system B. 

Lastly, we observe that as the number of particles in A grows, the current of the sub- 
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FIG. 3: (Color online) (a) Pearson coefficient vs. N. Squares and circles correspond to the scenario 
of one and two particles in B, respectively. Empty and filled symbols correspond to the t s and 77 
time scales, respectively. Inset shows the evolution of the current of two particles in B in units 
of C/h for N = 11. (b) Normalized Schmidt number 1C vs. time, (c) J vs. time. In (b) and (c) 
dashed-line: C/\U\ = 10 -4 , solid line: C/\U\ = 1/10 .The other parameters are the same as in 
FigJU 



system with few particles tends to evolve towards an equilibrium value very close to zero, 
corresponding to J ~ 1 (see Fig. EJ^c)). To understand this let us analyze the current of 
particles in the subsystem B 

J B = 7?TT^\ 5^ (m'\p B \m)(m\j B \m') (15) 
Tr(p B ) 

m,m' 

where \m) = &JJ0). Now, since the only non-zero matrix elements of the current are 
(m\ J B \m ± 1) then only off-diagonal elements of p B would contribute to the current. Nev- 
ertheless, it is known that in the limit ua 3> n B entanglement between particles in B and 



A results in a thermalization process; 
density matrix are reduced to zero 



19], where the non-diagonal elements of the reduced 
20|. Thus, the observation that J B equilibrates to zero 
is a direct consequence of the nature of the relaxation dynamics. Remarkably, the current 
could be used not only as an indicator of entanglement generation but also of how close the 
system is to thermal equilibrium. 

To conclude, we have studied the physical detection of entanglement for two sets of 
different atomic species moving in a ring-shaped optical lattice. The setup allows high 
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entanglement generation and also permits the estimation of the degree of entanglement by 
measuring the flux of particles. We have found that the generation of entanglement is highly 
correlated to the modification of the current in one of the subsystems. This current could 
then be used as a "witness" of entanglement generation in future experiments. Furthermore, 
it was shown that the same current can be used as a tool for the detection of thermalization 
between species of atoms. 
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also supported by FONDECYT project no 11110537 and MO is supported by FONDECYT 
project no 1100039. 
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